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We numerically study protocols consisting of repeated applications of two qubit gates used for 
generating random pure states. A necessary number of steps needed in order to generate states 
displaying bipartite entanglement typical of random states is considered. Numerics indicates that 
for a generic two qubit entangling gate the decay of purity is exponential with the decay time scaling 
as ~ n, implying that of order ~ n 2 steps are needed to reach random bipartite entanglement. 
We also numerically identify the optimal two qubit gate for which the convergence is the fastest. 
Perhaps surprisingly, applying the same good two qubit gate in addition to a random single qubit 
rotations at each step leads to a faster generation of entanglement than applying a random two 
qubit transformation at each step. 
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I. INTRODUCTION 

Entanglement is one of the resources which can be used 
to perform tasks not possible by classical means. Con- 
siderable effort has been put into understanding entan- 
glement properties of various states. While entanglement 
of certain classes of states, for instance of random states, 
is well understood, quantifying entanglement of a gen- 
eral quantum state is a rather difficult task, for a review 
see, e.g. In the present paper we are going to study 
optimal protocols for producing random states from ini- 
tial product states. By random states we mean pure 
quantum states which are eigenvectors of random unitary 
matrices distributed according to the unitarily invariant 
Haar measure. Random states have been the subject of 
considerable interest in the past. Entropy and purity of 
a subsystem have been studied in 0, y|, 0, [E B, 0, Hi • 
Distribution of the inverse of purity has been consid- 
ered in Q, the distribution of G concurrence (geometric 
mean of Schmidt coefficients)has been derived in [Ic| |. 
the distribution of purity i n 1 1 1| and the average values 
of Schmidt coefficients in [12|]. It has been shown that 
random states reproduce certain statistical properties of 
eigenstates of quantum systems whose classical counter- 
parts are chaotic very well. Even though a particular 
chaotic system has a well defined Hamiltonian a statis- 
tical description with random Hamiltonian, the so-called 
random matrix theory [13j |. has been very successful. 

One can ask if random states are also relevant for 
quantum information theory? Note that if description 
in terms of random states is applicable this usually sim- 
plifies the analysis [l4j |. Random states are a directly 
needed resource in certain procedures like quantum dense 
coding [l5| or remote state preparation [lri Il7| . One 
area at which quantum physics can do better than clas- 
sical is at computation, for instance, simulating quan- 
tum systems. If a state during quantum computation 
is not sufficiently entangled, efficient classical simulation 
is possible [H, [l9j]. Becouse random states are almost 
maximally entangled it is reasonable to expect that such 
random states will naturally occur during sufficiently 



complex quantum computation. In fact, we know that 
the evolution with chaotic systems will produce states 
whose statistical properties are well reproduced by ran- 
dom states. Therefore, quantum simulation of chaotic 
system is a likely candidate for problems which can be 
efficiently solved by quantum computer but not by clas- 
sical. For examples of such algorithms see, e.g. (20. |2l|. 
Note that sufficient entanglement is necessary but not 
sufficient for a quantum speed-up. An example are for 
instance quantum integrable systems which can also ef- 
ficiently produce a lot of entanglement [22], |23[ however, 
efficient classical simulation seems to be possible [24j |. 

An important practical question is how to produce ran- 
dom states efficiently? By efficiently we mean in a num- 
ber of steps that grows only polynomially with the num- 
ber of qubits n. A way to achieve this is already suggested 
by chaotic dynamics. Evolution with chaotic dynamics 
will produce states which are generic, i.e., random, from 
almost any initial state. It is well established that bipar- 
tite entanglement as measured by purity or von Neumann 
entropy increases linearly with time for chaotic systems 
and that the asymptotic saturation is reached after time 
which grows linearly with the logarithm of the Hilbcrt 
space dimension |25| . Typically, such chaotic Hamil- 
tonian can be written as a sum of n two qubit terms, 
and becouse entanglement saturates after time oc n, this 
means that of order ~ n 2 two qubit gates are needed 
to generate a random state - at least as far as bipar- 
tite entanglement is concerned. Our results will confirm 
these expectations. One should be aware that to pro- 
duce an arbitrary unitary transformation, and therefore 
a truly random state, an exponential number of two qubit 
gates is needed in general. However, if our criteria is just 
to reproduce bipartite entanglement of a typical random 
state, which is the case in the present paper, only poly- 
nomial number of gates is needed. The procedure to gen- 
erate random states is therefore the following: generate 
some pseudo-random sequence of gates, thereby produc- 
ing the so-called pseudo unitary operator, applying it to 
an arbitrary separable initial state. After sufficient num- 
ber of steps, i.e., applications of random gates, we will 
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end up in a random state. Such random protocol has 
been numerically studied in [26l [27j while a convergence 
to uniform Haar measure has been discussed in [28[ . The 
amount of interference produced by such protocols has 
been considered in [29]. Recently the question of how 
many gates do we need to obtain convergence has been 
attacked by analytical tools of Markovian chains [3(3] and 
a lower bound on the number of steps needed has been 
proved. In the present paper we are going to numerically 
calculate the exact convergence times for small chains 
as well as find a two qubit gate which will result in the 
fastest generation of bipartite entanglement. 



II. PROTOCOL 

Let us denote a state at time t by \ip(t)). The pro- 
tocol for generation of random entanglement consists of 
application of two qubit transformation at each time step 

m+\)) = u l3 {t)\m), (i) 

with Uij (t) being a two qubit gate acting on i-th and j-th 
qubits out of total n qubits. At each time step the pair 
of qubits i and j on which a gate acts will be drawn inde- 
pendently. We will consider three cases (couplings): (i) 
i and j are chosen randomly, that is gate can act on an 
arbitrarily separated qubits; (ii) i and j are neighboring 
qubits, that is a par (i, i + 1) or (i + and we take 
periodic boundary conditions, meaning that the first and 
the last qubits can also be coupled. Such coupling will 
be abbreviated nnPBC. (iii) similarly as in case (ii) we 
allow only nearest neighbor gates but with open bound- 
ary conditions, that is we do not allow a gate between the 
first and the last qubit. This will be abbreviated nnOBC. 
For the unitary matrix Uij we choose it to be a product 
of two independent single qubit unitaries Vi and Vj and 
two qubit gate Wij , 

U ij (t)=V i (t)V j (t)W ij , (2) 

Two qubit gate will be the same for all steps whereas sin- 
gle qubit unitaries Vi and Vj will be chosen according to 
the U(2) invariant Haar measure at each time step and for 
each qubit independently, that is they are from CUE(2) 
ensemble. Motivation for such a protocol is that from the 
experimental perspective two qubit gates are difficult to 
make and therefore we always apply the same two qubit 
gate, whereas single qubit transformations are relatively 
simpler and can be changed at each step. Without sac- 
rificing generality we will always choose the initial state 
to be a separable IV'(O)) = |00 . . . 0). All the statements 
about the convergence times in the paper thus pertain 
only to separable initial states. Choosing an entangled 
initial state can presumably lead to a faster convergence. 

The goal of entanglement generating protocol is to pro- 
duce states whose entanglement is as close as possible to 
that of random pure states, i.e., states drawn according 



to the invariant Haar measure on n qubits. Bipartite en- 
tanglement of a pure state is completely determined by 
its Schmidt coefficients pi, 

\ip) = ^^ilA® \<Pi)B, (3) 

i 

where \ipi)A and \<Pi)B are orthogonal and we assume pi 
are listed in nonincreasing order. Stating all Schmidt 
coefficients completely characterizes bipartite entangle- 
ment 31]. As a measure of closeness of our state \ip{t)) 
to a random state, we could for instance compare average 
values of Schmidt coefficients. For random pure states 
the average value of i-th largest Schmidt coefficient \ii 
has been calculated in In the case of a symmetric 
bipartite cut to first n/2 and last n/2 qubits it is given 
in an implicit form by 

2 cos (<fii) 

Ml = y/N ' 

ii±i^ = ^-i s in(2^), (4) 

where N = 2™/ 2 is the dimension of subspaces. For an- 
alytical treatment though, Schmidt coefficients of \i/j(t)) 
are not the simplest quantity to calculate. Therefore, in 
most of the paper we will rather use purity as a measure 
of entanglement. Purity I(t) is simply a sum of pf, or in 
terms of reduced density matrix 

m = tv A [ P \(t)] , PA (t) = tv B [\m)(m\] • (5) 

Here subscript A and B denote subspaces of first and the 
last n/2 qubits, respectively. For random states purity 
has been calculated @, 0] and is I(oo) x 2/N . In the last 
part of the paper we will also briefly mention results for 
the von Neumann entropy, 

S(t) = -tT A [p A (t)log 2 p A {t)]. (6) 

Value of S for random states is known 0, 0] and is 
S(oo) X n/2 — l/log c 4. Whenever we speak about 
purity or von Neumann entropy we will have in mind 
their approach to asymptotic values, i.e., I(t) — I(oo) 
and 5(oo) - S(t). 

The goal of the paper is to analyze how fast the purity 
decays to its asymptotic random state value, in partic- 
ular, how the decay time scales with n. In addition we 
are going to find a two qubit gate W © for which the 
(average) purity I(t) will decay the fastest. As we will 
see, this also means that the convergence of pi or S(t) 
is the fastest. In principle, to find the optimal gate, we 
would have to check the whole 15 parameter space of 
two qubit gates. However, each two qubit gate W can be 
decomposed as [H], [33[ 

Wij =Ai® Bj w(a x ,a y , a z ) A' t ® B' v (7) 
w(a Xl a y , a z ) = exp (i— [a x <7?a* + a, y afcr l j + a z afaj] j , 
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where A,A',B and B' are single qubit unitaries and 
a x,y,z are standard Pauli matrices. Two qubit transfor- 
mation w(a x , a y ,a z ) which is parameterized by three pa- 
rameters is called a canonical form of W. In our protocol 
([2]) single qubit gates V are random unitaries therefore, 
as far as entanglement is concerned, we only have to con- 
sider two qubit gates in its canonical form. Instead of a 
15 parameter space we only have to find the optimal gate 
among a 3 parameter set w. Furthermore, w(a x ,a y ,a z ) 
has certain symmetries. For instance, we have a relation 
—iaf<jjw(a x ,a y ,a z ) — w(a x + 2, a y ,a z ), and similarly 
for other a's. Rotating w by 7r/2 around y-axis, R y = 
cxp (— ier a 7r/4), we get R y R y w(a x ,a y ,a z )Rl <8> iJJ = 
w(a z , a y ,a x ). There is also a symmetry between 1 + a x 
and 1 - a x , iafw*(l + a x , a y , a z )a x } = w(l - a x , a y ,a z ), 
as well as between positive and negative parameters, for 
instance, afw(a x ,a y ,a z )a^ = w* (a x ,a y ,—a z ). Becouse 
entanglement produced by w* and w* is the same as by 
w, all these symmetries mean that it is enough to con- 
sider a's in the following range: 



1 > & x > 0>y > a z 



> 0. 



(8) 



When doing numerics we will pay special attention to 
three choices of two qubit matrices. For gate W we are 
going to consider a CNOT and an XY gate with the cor- 
responding unitary matrices in standard computational 
basis equal to 
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(9) 

form of CNOT has parameters 
whereas XY has canonical form 
u>xy(1, 1, 0). Note that the number of nonzero a's 
directly gives the number of CNOT gates needed 
to make such a gate out of single qubit gates and 
CNOT'spsi]. In fact, XY gate is equivalent (has the 
same canonical form) as the product of two CNOTs, 
w CNOT W -CNOT ) also called a DCNOT gate. Third case 

will be a two qubit uni tary Uij(t) chosen from U(4) 
invariant Haar measure [33], shortly U(4) gate. Note 
that only in this case we choose an independent two 
qubit gate at each step, in all other cases two qubit gate 
is the same for all steps. This last choice of U(4) gate 
will serve for comparison. Naively, one might think that 
such random two qubit transformation at each step will 
produce entanglement in the fastest way. This is not 
so though. The optimal gate will turn out to be the 
XY gate. Analytical calculation of purity for general 
gate w(a x , a y ,a z ) and our protocol @ is rather difficult. 
Things simplify though for gates W which preserve the 
Pauli group. 



III. SPECIAL CASE : MARKOV CHAIN 

Let us expand a pure state density matrix over prod- 
ucts of Pauli matrices, 



\m)(m\ 



(10) 



where cr" 1 denotes Pauli matrix Ui S {0, x, y, z} acting 
on i-th qubit, with the convention <r° = 1. We use a 
short notation a — (a±, . . . , a„). With the expansion 
(|10[) purity ([5]) is now simply given by 



1(f) 



1 

iV2 



a={a A B } 



(11) 



where the summation runs over all a which have identity 
on the subspace B, i.e., aj = 0, j = n/2 + 1, . . . , n. To 
obtain the decay of purity we therefore have to calculate 
the time dependence of c 2 (i). Averaging over U(2) in- 
variant Haar measure of single qubit matrices Vi, Vj @ 
we get the transformation law 

C a(t + 1) = VpjVyjRaiJiRajJjCpit)^^), (12) 

/3,7,<5 

where v a ,p is denned by Wij<r a W-j = J2p v a,i3< jl3 an d R 
is a 4 x 4 matrix obtained from averaging over random 
single qubit gates and is equal to 



R = 



/l 0\ 

± i i 
? ? ? 



vo I S ! 

V 3 3 3 



(13) 



If Wij preserves the Pauli group, i.e., if it transforms 
products of Pauli matrices into a product of some other 
Pauli matrices, that is if v^p . )7i7 . is nonzero and equal to 
±1 or ±i only if 7^ = [3[ and jj = (3j, the transformation 
([12")) can be simplified to 

ci(t + l)=J2(Mg ) ) a , c 2 fj (t), (14) 
(M^f) a ,i3 = R ait 0>.R aj> i3'.- 

(2) 

Markov matrix of dimension 16 x 16 involves only 

i-th and j-th qubits. To get the transformation on all 
n qubits we have to average over all couplings, resulting 
in Hi 



c 2 (i+ 1) = Mc 2 (t), 



M 



r(2) 

L ^ " y 

ij 



(15) 



if L is the number of all couplings, i.e., number of al- 
lowed pairs of qubits i and j. As mentioned, we will con- 
sider three different couplings, random i and j, nnPBC 
and nnOBC. Let us illustrate how f3' are determined 
from (3 (for instance, in ea. (|15|) ) on the example of a 
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CNOT gate. All 16 different products of Pauli ma- 

trices <7j ct,- 3 can be numbered by i = pj + 4 • pV 
Transformation of all 16 products can now be stated 
by giving the transformations of x's. For instance, the 
transformation Ucnotc x ^^cnot = aV<jZ can be sim- 
ply stated by x — 6 going to x' — /?'• + 4 • = 
11. All 16 x' for CNOT can be written in a vector 
x' = (0,1,14,15,5,4,11,10,9,8,7,6,12,13,2,3), denot- 
ing the transformation x — (0,1,2,..., 14, 15) — > a/. It 
turns out that XY gate also preserves the Pauli group. 
Transformed a;'s for XY gate can be written in a vec- 
tor (0, 11, 7, 12, 14, 5, 9, 2, 13, 6, 10, 1, 3, 8, 4, 15). Marko- 
vian description (1 1 5[) greatly simplifies the analysis and 
this formulation [30( will be used to calculate the decay 
of purity. Purity is given by a sum of certain c 2 ([TT]) 
and therefore its asymptotic decay will be determined by 
the second largest eigenvalue of the matrix M. It can be 
shown that M (fT5"|) has two eigenvalues equal to 1. One 
corresponds to the identity operator on all qubits, and 
the other to the ergodic density uniform on the rest of the 
space. A third eigenvalue, A3 = 1 — A, then determines 
the asymptotic decay of purity as I(t) — I (00) x (1 — A)'. 

In addition to two qubit gates that preserve the Pauli 
group, Markov dynamics of the form c 2 (i + 1) = Mc 2 (t) 
can also be written for the case when at each step we 
choose an independent U(4) gate. Averaging over U(4) 
group gives M' 2 ' which preserves identity, ofcr° — > of cr°, 
while it uniformly mixes all other 15 possible products 
oT a 7 ■ Matrix elements are therefore [M^ (UA)] ,o = 1 
and [MW(U4)] X , X , = 1/15 if x,x' € {1,...,15}. 

Most of two qubit gates w(a x ,a y ,a z ) though do not 
preserve the Pauli group and therefore transformation of 
c 2 (t)'s can not be written as a Markov process (fT"5|) . In 
such cases we will have to use a direct numerical sim- 
ulation of our protocol, averaging over many different 
realizations to obtain the decay of I(t). 

In the next section we are going to present numerical 
results on the convergence rates. Most of the time we are 
going to focus on purity. At the end of the section we will 
see that other quantities like the von Neumann entropy 
or the Schmidt coefficients give essentially the same in- 
formation. Convergence rate will be obtained by several 
different methods. When looking for the optimal gate, a 
direct simulation of the protocol is used to calculate I(t) 
at some fixed time for different gates, thereby identify- 
ing the optimal one. Central quantity of interest is the 
scaling of the convergence rate with n. In this respect 
we numerically calculate the gap of the corresponding 
Markov chain for different numbers of qubits n. We find 
that the gap in all cases scales as A ~ 1/n. Finally, we 
use the obtained gaps A to predict the asymptotic decay 
of purity. 



IV. NUMERICAL RESULTS 
A. Optimal gate 

First we want to find an optimal two qubit gate W @ 
which will result in the fastest possible decay of purity 
([5]). We have already discussed that due to U(2) invari- 
ance of single qubit gates Vi we can limit the study to 
two qubit gates in the canonical form J5]) with the pa- 
rameters a x ,y,z in the range 1 > a x > a y > a z > ([5]). 
Becouse Markov description for the transformation of c 2 
is not possible for general parameters we have to resort to 
direct numerical calculation of purity. For each set of pa- 
rameters we calculated the average purity I(T) at some 
fixed time T, from which we then deduced the expected 
decay rate k of purity, assuming the dependence 

I(t) - I(oo) = eacp (-«*/«)• (16) 

Optimal gate is then the one with the largest k. We 
always took n = 8 qubits and averaged over 1000 protocol 
realizations. Time T at which we calculated k was 30 for 
random i — j coupling and 50 for nnPBC and nnOBC 
couplings. Note that becouse in each case n and T were 
fixed we could instead of k simply use / as a criterion to 
determine the optimal gate. The reason to use k is, as 
we will see later, that such theoretical form of purity is 
predicted in most cases of Markovian dynamics. 




FIG. 1: (Color online) The dependence of the decay rate of 
purity — k (|16[) on three parameters a x , y ,z for a protocol with 
random i—j coupling. Cross-sections for four different a z are 
shown. Fastest decay (« ~ 1-2)) is obtained for gates of the 
form w(l,a y ,0). This includes XY gate at w(l, 1,0) as well 
as CNOT gate at w(l, 0, 0). Dashed triangle shows the set of 
parameters fulfilling Eq. Q. 

In Fig. Q]we show the results in the case of coupling be- 
tween random pairs of qubits. It seems that the fastest 
decay of purity is achieved for a continuous family of 
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two qubit gates of the form w(l,a y ,0) (and all vari- 
ants obtained by symmetries). This includes XY gate 
with w(l,l,0) as well as CNOT with w(l,0,0). Gate 
w(l, 1, 1), which is equivalent to a SWAP gate, does not 
produce any entanglement. 




FIG. 2: (Color online) The dependence of the decay rate of 
purity —k (|16[) on three parameters a X:VtZ , similarly as in 
Fig. [1] On the left is for nnPBC coupling and on the right for 
nnOBC. 

In Fig. [2] we show similar plots for nnPBC and nnOBC 
coupling. Apart from the numerical values of k the over- 
all dependence on a's looks similar to the one for random 
i—j coupling. There is one notable difference thou. Now 
the XY and CNOT gate do not generate entanglement 
equally fast. For nnPBC for instance, k « 0.7 for XY 
and 0.5 for CNOT. XY is clearly better for generation 
of random entanglement. Similar conclusion is reached 
for nnOBC coupling. Interestingly, comparing nnPBC 
and nnOBC cases we can see that k's for the later are 
by a factor w 2 smaller than for nnPBC. Remember that 
the only difference between nnPBC and nnOBC is in the 
coupling between the first and the last qubit, i.e., in the 
term Uo, n -i, which is absent in nnOBC. This sole differ- 
ence results in nnOBC being by a factor of w 2 slower 
in producing random entanglement. One can intuitively 
understand this factor in the following way. The aver- 
age distance between qubits in subspaces A and B is in 
the case of nnPBC (qubits on a circle) about two times 
smaller than for nnOBC (qubits on a line). Therefore, be- 
couse the decay rate for purity (fTo| scales inversely with 
n one can expect a factor of 2 between the two cases. 



B. Markov chain 

In previous section we have identified two interesting 
gates, XY and CNOT, which are optimal for certain cou- 
plings. As we have seen in section IIIII random proto- 
col can be described by Markov chain for both of these 
gates. This has several advantages. In comparison to 
numerical calculation of purity averaging over protocol 
is in Markovian formulation exact. Provided the sec- 
ond largest eigenvalue (i.e., the largest smaller than 1) of 
Markov matrix M is nondegenerate, 1 — A, the decay of 



purity will be for small gaps A given by 

7(t)-/(oo) = exp(-t/r), r = 1. (17) 

Therefore, knowing the gap we will know the decay rate. 
Using analytical techniques one can actually bound the 
gap. This has been done in [30l | where they proved that 
the gap is A > g^j^irn f° r CNOT gate and random i — j 
coupling. As the exact analytical calculation of the gap 
seems too difficult, even for a relatively simple Markov 
chain with U(4) gate, we are going to numerically cal- 
culate the values of the gap for different n and differ- 
ent couplings. In all cases we are going to consider XY, 
CNOT, and U(4) gate. Disadvantage of Markovian de- 
scription is that the matrix M is of rather large size 4", 
i.e., their dimension equals to the square of the Hilbert 
space dimension of pure states. 

Results of numerical calculation of the gap for a chain 
with random i — j coupling are in Fig. [3] We use the 
Lanczos method [35[ to calculate few largest eigenvalues 
of the Markovian matrix. For n < 13 we use the original 
Markov chain M of dimension 4™ while for lar ger n a 
reduced chain is used whose size is only 2™, see [3(| for 
details. Becouse the size of the matrix grows exponen- 
tially with n we are limited to relatively small number 
of qubits. In accordance with Fig. [T]the gap is the same 
for XY and CNOT gates (same within numerical preci- 
sion). Fitting 1/n dependence to the values of the gap 
we get A(XY,ran-ij) = 1.47/(n + 2.15) for XY and 
CNOT gate and A(U4,ran-ij) = 1.33 /(n + 2.50) for 
U(4) gate [38|. Perhaps surprisingly, the U(4) is about 
10% slower than XY gate. Doing random two qubit gates 
at each step is therefore not as efficient in producing ran- 
dom entanglement as doing random single qubit gates 
and a fixed good two qubit gate (XY or CNOT)! Nu- 
merical result for the gap in the case of a CNOT gate 
improves the analytical lower bound A > 4/9n(n — 1) 
proved in 30]. For nnPBC results are shown in Fig. [4] 
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FIG. 3: Dependence of the inverse gap, 1/A, on n for Marko- 
vian chain with random i — j coupling and XY (triangles), 
CNOT (circles) and U(4) (squares) gate. Full lines are best 
fitting linear function, see text for details. 
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Here full symbols correspond to numerical calculation of 
the gap on a 4™ dimensional Markov chain while empty 
symbols are gaps indirectly determined from the purity 
decay (|T71) . The gap now scales as A(XY,nnPBC) = 
0.45/(n- 2.50), A(U4,nnPBC) = 0.36/ (n- 2.67) and 
A(CNOT,nnPBC) = 0.28/(n - 3.01) [39]. Here we 
see that XY is the optimal gate, while U(4) is better 
than CNOT gate, being the worst of the three. For 
nnOBC very similar results are obtained as for nnPBC 
and we will only list the fitted dependence of the gap, 
A(XY,nnOBC) = 0.23/(n - 3.01), A(U4, nnOBC) = 
0.19/(n - 3.12), A(CNOT,nnOBC) = 0.15/(rc - 2.98). 



50 
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n 

FIG. 4: Dependence of the inverse gap 1/A on n for 
nnPBC coupling and XY (triangles), CNOT (circles) and 
U(4) (squares) gates. Full lines are best fitting linear func- 
tions, see text for details. 



C. Purity decay 

While the asymptotic decay of purity is guaranteed to 
be I(t) — I(oo) >c exp (— iA) if the eigenvalue 1 — A of 
Markovian matrix is nondegenerate, the decay can be 
more complicated in the case of degeneracies. It is found 
that for many Markov chains one has the so-called cut- 
off phenomenon [36[. Briefly, we say there is a cutoff in 
a Markov chain if the distance between the asymptotic 
crgodic distribution and the distribution after t steps su- 
denly drops from 1 to at some cutoff time. The sharp- 
ness of this transition increases in the limit of large state 
space size (in our case for n — > oo), for exact definition 
see [36| . Although precise mathematical conditions lead- 
ing to the cutoff are not known [361 ], it generally occurs 
due to multiplicity of the largest nontrivial eigenvalue 
with the degeneracy increasing with increasing n. Let 
us illustrate by an example: suppose we have n times 
degenerate largest (non 1) eigenvalue equal to 1 — 1/n. 
Purity will then decay as I ~ nexp(— t/n). From this 
we see that / will be equal to exp (— c), where c is some 
fixed number, at time t = n(lnn + c). However large 
c we choose, in the limit n — > 00 purity will be small 



/ ~ exp (— c) at time nlnn. In other words, for times 
slightly smaller than n In n purity will be very large while 
it will be very small for times slightly larger. There is a 
sudden cutoff at time nlnn. On the other hand, if the 
eigenvalue is nondegenerate, purity is going to decay as 
/ ~ exp (—t/n). It will reach small value exp (— c) at time 
nc. This time now increases with increasing c. There is 
no cutoff. If we want / to be smaller we have to increase 
t. One can also look at time derivative of purity. In 
the case of a cutoff the absolute value of the derivative 
(steepness of the curve I(t)) diverges at some fixed value 
of /, while it goes to zero if there is no cutoff. 
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FIG. 5: Decay of purity for n = 16 and random i—j coupling 
and U(4) and CNOT gates, dashed curves. Theoretical de- 
cay (|18(l is shown with the full curve. Asymptotically purity 
decays as exp (—tA) (|17|) with the gap predicted from Fig. [3] 
For smaller times thou degeneracies are important, reflected 
in a steeper decay of purity (|18|l . In the limit n — > oo the 
cutoff is predicted. 

To identify a possible cutoff one therefore has to look 
at multiplicity of the largest eigenvalue. In all cases we 
studied in Sec. IIV 51 the largest eigenvalue 1 — A is non- 
degenerate. However, large degeneracy of an eigenvalue 
that is very close to 1 — A could also possibly cause a 
cutoff. We have therefore numerically checked for de- 
generacies of the first three largest eigenvalues (here we 
mean first three nontrivial; we do not count a doubly de- 
generate eigenvalue 1) for gates XY, CNOT, and U(4). 
In all three cases we found degeneracies only for random 
i — j coupling, while there were no degenerate eigenval- 
ues for nnPBC and nnOBC couplings. For XY and U(4) 
gate and random i — j coupling a 3rd eigenvalues is n — 1 
times degenerate, while for CNOT gate a 2nd eigenvalues 
is n — 1 times degenerate. This leads us to speculate that 
there is a cutoff for all three gates in the case of a random 
i — j coupling [30] while there is probably no cutoff for 
nnPBC or nnOBC coupling. 

To check the above prediction about the cutoff as well 
as to verify the asymptotic decay of purity, we have per- 
formed numerical calculation of purity decay for larger 
n, where the cutoff phenomenon should be visible. We 
simulated our random protocol, averaging over many re- 
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alizations, thereby obtaining I(t). In Fig. [5] we show the 
decay of purity for n = 16 qubits and random i — j cou- 
pling. As we have seen, the largest eigenvalue is in all 
cases nondegenerate while the 2nd or 3rd eigenvalue is 
n — 1 times degenerate for random i — j coupling. There- 
fore, for short times the decay of purity will be given by 
the largest as well as by the mentioned degenerate eigen- 
value. We fitted numerically obtained curves in Fig. [5] 
with the following dependence 



I(t) - 7(oo) 



(1 + a ■ cxp (-b ■ t/r)) 



Numerical results presented show that the decay time 
of purity asymptotically scales as r ~ n. This holds 
also for gates that can not be described by Markovian 
chain (data not shown). Becouse random state value of 
purity is ~ 1/Y, this implies that purity will decay to 
its asymptotic value I(t) — -T(oo) ~ 1/2™ after number of 
steps scaling as t ~ n 2 . 



D. Measures in addition to purity 



exp(-Vr), (18) 



where r = 1/A is given by the gap of the largest eigen- 
value (|17p and is determined from linear fitting lines in 
Fig. [3l while two fitting parameters a and b take care 
of the degeneracy and the gap of the largest degenerate 
eigenvalue, respectively. For the CNOT gate and n = 16 
we get a = 2.80 and b = 0.75, while for U(4) gate we get 
a = 2.57 and b = 0.75. Decay times r (fT8|) used were the 
ones predicted from Fig. [3] and are r = 12.3 for CNOT 
and t = 13.9 for U(4) gate. Similar fitting parameters 
as for CNOT gate are obtained also for XY gate (data 
not shown). Important point is that the parameter a, 
connected with the degeneracy, increases with n. For in- 
stance, for U(4) gate it is a = 0.78 for n — 12, a = 2.57 
for n = 16 and a = 6.4 for n = 18. This signals the 
emergence of a cutoff for random i — j coupling and large 
n. 
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FIG. 6: Same as Fig. [5] but for nnPBC coupling. Purity is 
well described by theoretical I(t) — I(oo) = exp(— t/r) (|17l) . 
with r predicted in Fig. [4] (i.e., no fitting), in the whole range 
of times. No cutoff is expected. 

On the other hand, things are quite different for 
nnPBC or nnOBC coupling. Results of numerical simu- 
lation for nnPBC are in Fig. [6] (for nnOBC data are very 
similar). As one can see, purity is well described by a sin- 
gle exponential function ([17) with the predicted t = 29.9 
for XY, t = 37.0 for U(4) and t = 45.9 for CNOT gate. 
Therefore, according to numerical results for sizes upto 
n = 20, we can predict that there is probably no cutoff 
phenomenon for nnPBC and nnOBC couplings. 
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FIG. 7: Similarly as I(t), von Neumann entropy ([6} also de- 
cays on the same time scale as purity. All is for nnPBC, gates 
XY, CNOT and U(4) and n = 16. 

In the core of the paper, due to its simplicity, we have 
used purity I(t) as a measure of entanglement. Here we 
are going to show that other quantities, von Neumann 
entropy and Schmidt coefficients, for instance, also decay 
on the same time scale r = 1/A as purity. The only 
difference is that the functional dependence of the decay 
is more complicated than a simple exponential function. 
First we show the results for the von Neumann entropy 
S(t) ©. In order to show the scaling with decay time 
t we plot in Fig. [7] the dependence of S on t/r. As we 
can see, decay of S(t) is indeed described by a universal 
form S(oo) — S(t) = f(t/r), with the function / visible 
in figure. It slightly depends on the type of the coupling, 
i.e., whether we have a random i — j, nnPBC or nnOBC, 
but is independent of the gate. 

To fully specify bipartite entanglement a single quan- 
tity, like purity or von Neumann entropy, is not enough. 
One has to specify all Schmidt coefficients /i,-. As a fi- 
nal test we checked how individual Schmidt coefficients 
converge to those of random pure states. Theoretical 
prediction for random pure states has been calculated 
in [12| and is given in an implicit form by Eq. ^ . Again, 
Fig. [5] shows the convergence on the same time scale as 
that of purity. Functional dependence looks Gaussian for 
nnPBC coupling while it is more complicated for random 
i — j coupling, possibly due to degeneracies. In Fig. [5] we 
show the convergence of individual /z?. 
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FIG. 8: Schmidt coefficients also converge on the same time 
scale r as purity. We show how £\ (t) — fii(oo) \ depends 
on t/r. On the left is for random i — j, on the right for 
nnPBC coupling (full line is a Gaussian fit). In both plots 
three curves almost overlapping are for CNOT, XY, and U(4) 
gates, n = 16. 
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FIG. 9: Convergence of squares of Schmidt coefficients Hi(t) 
(i.e., eigenvalues of the reduced density matrix) for n — 16 
and XY gate with nnPBC. Theoretical decay time is r = 29.9. 
Full thick line is theoretical prediction ^j?(oo) ©. 



V. DISCUSSION 



We have numerically studied protocols for generation 
of generic entanglement as represented by random pure 
quantum states. At each step of the protocol a fixed 
two qubit gate and two independent random single qubit 
unitaries are applied. We have calculated the decay rate 
of purity which in all cases studied grows linearly with 
the number of qubits n. For certain two-qubit gates 
Markovian description is possible. Convergence rate is 
in such cases determined by the size of the gap which is 
numerically found to scale as 1/n, improving the analyt- 
ical bound in [3(J. An optimal two qubit gate is identi- 
fied which produces random entanglement in the small- 
est number of steps. This optimal gate is for all different 
couplings considered XY gate, also called DCNOT gate, 
generated by the Heisenberg XX interaction. Depending 
on the coupling XY gate can be as much as 60% faster 
than CNOT gate. Interestingly, applying a random two 
qubit gate at each time step is slower than applying a 
fixed good two qubit gate in addition to random single 
qubit unitaries. For coupling between random qubits we 
predict the cutoff phenomenon, while there is probably 
no cutoff for the protocol where we apply two qubit gates 
only to nearest neighbors. 
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